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Abstract: An outline is given of an extended perturbative solution of Euclidean QCD which system- 
atically accounts for a class of nonperturbative effects, while still allowing renormalization by the 
perturbative counterterms. Euclidean proper vertices T are approximated by a double sequence 
where r denotes the degree of rational approximation with respect to the spontaneous mass scale 
Aqcd, nonanalytic in the coupling g 2 , while p represents the order of perturbative corrections in g 2 
calculated from r[ r ' l - rather than from the perturbative Feynman rules r(°) pert - as a starting point. 
The mechanism allowing the nonperturbative terms to reproduce themselves in the Dyson-Schwinger 
equations preserves perturbative renormalizability and is intimately tied to the divergence structure of 
the theory. As a result, it restricts the self-consistency problem for the r[' r,0 l rigorously - i.e. without 
decoupling approximations - to the seven superficially divergent vertices. An interesting aspect of the 
solution is that rational-function sequences for the QCD propagators contain subsequences describing 
short-lived elementary excitations. The method is calculational, in that it allows the known techniques 
of loop computation to be used while dealing with integrands of truly nonperturbative content. 



1 Generalities and Notation 



1.1 Nonperturbative Quantities 

One of the more important insights to have emerged from two decades of study of the large-order 
behavior and summability of perturbation expansions [1] has been that for renormalizable but not 
superrenormalizable field theories (typically in four dimensions), the perturbative series of correlation 
functions T(g 2 ) around g 2 = is fundamentally incomplete, in the sense that it does not allow unique 
reconstruction of those functions even in principle. Perturbation series, in most cases of interest, are 
divergent asymptotic series which at fixed positive g 2 allow a function to be estimated with a finite 
accuracy. A typical error estimate is 
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as can be inferred from the behavior of the perturbative coefficients rto pert generally found in quan- 
tum field theories [1]. For superrenormalizable theories, typically in dimensions D < 3, analyticity 
properties of the T's in the complex g 2 plane can be established outside of perturbation theory that 
are strong enough to conclude that a function sufficiently analytic and satisfying this bound must in 
fact vanish, and that T's can be reconstructed uniquely from their asymptotic series by resummation 
techniques. By contrast, for renormalizable theories in D = 4, the necessary analyticity is positively 
known to be violated [2], and T's may therefore differ from their perturbative series by terms ex- 
ponentially small near g 2 = that can "duck under" the bound (1.1) and remain invisible in the 
expansion. 

The general mathematics of semi-convergent series does not say more about the missing terms, 
if any, that have no expansion in g 2 in even a formal sense : one must turn to physics for clues. 
Asymptotically free theories are special in that there one has a-priori knowledge, from renormalization- 
group (RG) analysis, of the presence of one important class of terms just allowed by the bound (1.1). 
These are the terms involving the RG-invariant spontaneous mass scale [3], 



A 2 (g 2 (y),^ R = u 2 exp 
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(1.2) 



where v is the arbitrary mass scale introduced by renormalization in a scheme R. With asymptotic 
freedom, 

A 2 = ^exp(-— T i— -[l + 0(g 2 )}} , (1.3) 



(where /3b, defined by eq. (1.20) below, is positive), the archetypal non-analytic g 2 dependence of the 
"exponentially small" form. It hardly needs emphasis that the term "exponentially small" refers to 



1 



the hypothetical behavior of such a quantity in a purely formal, unphysical limit - the limit where one 
sends g 2 — > without a corresponding increase in v in the sense of RG flow. In the commonly used 
(dimensional) schemes R, A# of course is not small, and in fact sets the scale of hadronic masses. 

It is also clear from the context that when speaking of A 2 -dependent terms we are not alluding here 
to the relatively weak A 2 dependence obtained by the standard RG process of introducing a running 
coupling constant a(k 2 /A 2 ), with asymptotic expansion 

5(2) - 1 <i+o(±)\ (,=£»/). a,. 
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since this arises purely within the perturbation series by leading-logarithms (LL) resummation. By 
contrast, consider some quantities presumed to be nonperturbative in the deeper sense mentioned 
above, such as the vacuum expectation of the trace anomaly [4] of the energy-momentum tensor in 
massless QCD. Since this is known to be an RG invariant and of mass dimension four, it must be of 
the form 



< 



^ G» v G£" 



2q „ r,.V ; . (1.5) 

where ca is a pure number. As another example, again in massless QCD, consider dynamical forma- 
tion of a color-singlet glueball state, whose mass defines an RG-invariant scale. The connected and 
amputated four-gluon amplitude takes the form, 

T abcd (pi Pa) = S ab S cd {pi+p2)2 + CBA2 (1-6) 

+ (crossed terms) + (regular terms) , 

with eg again a pure number. Yet another example is furnished by instanton effects. There, n- 
instanton contributions to physical quantities typically appear with factors exp(— nS\/ g 2 ), where the 
instanton action S\ is in general not an integer multiple of the (4tt) 2 / (3q of eq. (1.3), so that fractional 
powers of A seem to be present. However, strong arguments have been given by Miinster [5] to the 
effect that this is really an artefact of the dilute-gas approximation used in most instanton calculations, 
and that by taking account of the denseness of the instanton system one is led back to integer powers 
of A in the properly summed multi-instanton series for physical quantities. These examples indicate, 
as does experience with operator-product expansions, that the A dependence of truly nonperturbative 
quantities will typically be of a polynomial or rational form. It is this type of dependence which can be 
clearly separated, without double-counting problems, from the "resummed-perturbative" one of (1.4), 
and for which g 2 and A 2 can be handled formally like two different parameters, (although in the end 
one will always seek to use LL resummation to eliminate g 2 and deal only with a single, .KG-invariant 
parameter, A 2 ). 
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The present paper outlines [6] , for an Euclidean and asymptotically free gauge theory exemplified 
by QCD, a systematic extension of perturbation theory designed to account for nonperturbative effects 
as reflected in a rational dependence, or more generally a dependence representable by sequences of 
rational functions, of correlation functions on A. Since the approximating sequence is to provide a 
genuine extension, rather than a resummation, of the perturbative one, it is natural that it should 
take the form of a double sequence, T^ r ' p \ as described in sect. 2 for the simplest case of a one- 
component vertex function depending on a single squared Euclidean momentum (the two-point vertex 
of transverse gluons); sect. 3 and appendix A describe the relatively straightforward extension to the 
more complex (three- and four-point) basic vertex functions. In this double sequence, the integer r 
will characterize a certain degree of approximation with respect to A, while the "perturbative" index 
p still counts explicit powers of g 2 (y). That r cannot simply be an index counting powers of A has 
to do with the fundamental difference of dimensionality: with A being a mass, as opposed to the 
dimensionless g 2 , an expansion in powers of A would (for fixed mass dimension dr of T) inevitably be 
a large-Q 2 expansion in powers of A 2 /Q 2 , where Q 2 is a typical squared Euclidean four-momentum of 
r, and such an expansion would never provide an adequate representation of T in the region Q 2 < A 2 , 
whereas the dynamical equations determining T's contain loop integrals whose evaluation requires a 
systematic approximation of those T's over the entire momentum range. Of necessity, the index r 
will therefore have to refer to a sequence of more global approximants, whereas for g 2 , which enters 
the dynamical equations only parametrically, a local approximation around g 2 = is possible. The 
peculiar asymmetry between the two "directions" of the sequence is therefore not an arbitrary choice, 
but is rooted in the very nature of the spontaneous scale (1.2). 

On the other side, the meaning of p will also be subtly different from what it is in the purely 
perturbative context, since for p > 1 it will refer to corrections calculated from diagrammatic building 
blocks rl r '°], rather than from the standard Feynman-rules vertices Y^ pert implied in eq. (1-1). 
In other words, as compared to perturbation theory, the approximation will be able to "correct its 
own zeroth order". Determination of the set {r^' !}, the nonperturbatively modified vertices of zeroth 
perturbative order, will then lead to the self-consistency problem outlined in sect. 4, which is one of the 
characteristic new features of the scheme. We discuss the remarkable fact that this self-consistency 
problem becomes rigorously restricted, not by any "decoupling" approximations but by the very 
nature of a mechanism tied to the divergence structure of the theory, to a small finite set of vertices, 
and we indicate how perturbative renormalization continues to function in the framework of this 
mechanism. Finally, section 5 has a number of comments on the special description emerging in 
this context of the elementary QCD excitations. Apart from a few glimpses provided by the sample 
calculation of sect. 4, this outline does not dwell at all on either the bulk of the actual loop and self- 
consistency computations, nor on the interesting, but nevertheless dynamically secondary, question of 
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approximate nonperturbative saturation of the Slavnov- Taylor identities, since both subjects are of 
an algebraic lengthiness that requires separate presentation. We do, however, comment occasionally 
on the relation with the work of refs. [7], which the present scheme allows to put in perspective 
as a lowest stage (although with somewhat arbitrary technical simplifications) of a more systematic 
sequential approximation. 

1.2 Notation and Conventions 

The following notation will be used. One considers correlation functions generated by the standard 
gauge-theory action in D = 4 — 2e Euclidean dimensions, 




where Cy and Cp are Lagrange densities, respectively, of a set of SU(Nc) gauge- vector (gluon) fields 
and of a set of minimally coupled fermion (quark) fields coming in Np flavors F = u,d,s..., while 
the term 

C G = ^[d^A^x)] 2 + c a (x){[-5 ab d^ + g f abc A^x)]d^}c b (x). (1.8) 

comprises standard covariant gauge-fixing and Fadde'ev-Popov (FP) terms. The coupling factor go = 
v^go, with vq another mass scale, is used to keep the bare gauge coupling go dimensionless. We will 
focus in particular on the connected, amputated and one-particle irreducible correlation functions 
(proper vertices) with N external lines in Euclidean momentum space, 

rV(W), {k} = {k 1 ,k 2 ,...k N \k 1 + k 2 + ... + k N = 0}. (1.9) 

Where necessary, the ./V external lines will be labelled more specifically by V for a vector (gluon) line 
- occasionally detailed further as T or L for a transverse or longitudinal gluon, respectively -, by G, G 
for ghost and antighost lines, and by F, F for quark and antiquark lines. Thus Ttlv will be a purely 
gluonic three-point vertex with one transverse, one longitudinal, and one generic gluon line, r 'pyyp 
a quark-gluon-gluon-antiquark four-point vertex, etc. The numbers of the three types of lines will be 
written ny,nc,nF, so that 

n v + n G + J2 n F = N - ( L1 °) 

F 

In particular, the set of proper two-point vertices 

E 2 := { T vv = -D~\ T GG = -F>-\ {^ = -5^ = 1...^}}, (1.11) 

consists of the negative inverses of the vector propagator D, ghost propagator D, and fermion propa- 
gators Sp, while the set 

T 3 := { T 3V , T GVG , {T FVF \F = 1 . . . N F } } (1.12) 
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comprises the basic three-point interaction vertices. The dynamics implied by the action (1.7 - 1.8) are 
then embodied in the Dyson- Schwinger (DS) equations [8], an infinite, hierarchical system of coupled 
integral equations for the functions (1.9) (see [9] for their specific form in QCD), which in a condensed 
notation take the form 

rv = r$ )pert + ^£^ 2 ^[r 2 ,r 3 ,..., t n ,t n+1 ,t n+2 ] . (1.13) 

The perturbative zeroth-order or bare vertices, r^ pcrt , are given by the standard Feynman rules for the 
action (1.7 - 1.8), while <&n denotes a set of nonlinear dressing functionals, defined by loop integrals 
over combinations of T's, and the notation emphasizes that each such loop integral is preceded by at 
least one power of the bare gauge coupling, g%. 

It is important to keep in mind that a fundamental dichotomy is introduced into the set of vertices 
(1.9) through the renormalizable divergence structure of a QCD-like theory. A small finite subset, 
consisting of the superficially divergent or basic vertices, 

Tsdiv := { T 2 , T 3 , T AV } , (1.14) 

are distinguished by the fact that loop integrals in their $jy functionals of eqs. (1.13) have a non- 
negative value of the overall degree of divergence in D = 4 [10], 



5n = 4 — ny — 




(1.15) 



so that each of these need their own specific renormalizations, which can be performed at least per- 
turbatively [11]. Related to this is the fact that it is precisely for T S( n v that the bare terms r(°) pert in 
eq. (1.13) are nonzero. By contrast, the remaining, infinite set of superficially convergent or higher 
vertices, 

^sconv := { ^GVVGi ^ FVVFi ^ FFFFi ^W, • • •} (1-16) 

have loop integrals with a negative <5^, and therefore exhibit no typical divergences of their own, 
but at most subdivergences representing corrections to the basic vertices (1.14): when rewritten as 
dressed- skeleton expansions, i.e. in terms of fully dressed and renormalized basic vertices (1.14), their 
loops are actually convergent. 

While testable consequences of the theory are mostly contained in (the color-singlet channels of) 
the higher Green's functions (1.16), these cannot be calculated in a truly systematic way without first 
studying, and renormalizing, the basic vertices (1.14): a disproportionate amount of theoretical effort 
must be directed towards a class of amplitudes that contain next to nothing in observable physics. 
(An exception may be the four-gluon vertex T^v in (1.14), which may develop glueball poles.) The 
approach described below is, at its present stage, concerned exclusively with the DS equations (1.13) 
for the seven superficially divergent vertices (1.14). 
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The usual perturbative solution to eqs. (1.13) is obtained by straightforward iteration around 
p(o)pcrt anc j applying, a t each step, a renormalization scheme R, which among other things elimi- 
nates go in favor of a renormalized coupling g{y) depending on the arbitrary scale u : 



-ppert _ r [p]pert, F [p]pert _ -p(0)pert ^ 

1 N - p 1 ™^ 1 N ' 1 N - 1 N +2^ 



p'=l 



4-7T 



-,(p')pert 
N 



(1.17) 



Here the radiative corrections r^ pcrt (p > 1) are computed iteratively from the zeroth-order solu- 
tion r^ pcrt , i. e. from the standard Feynman rules. For example, the first iteration is described 
schematically by 




p(0)pert 



R,v 



4tt 



rj )pcrt + o( 5 4 ) . 



(1.18) 



For R, we always have in mind a dimensional-regularization-plus-minimal-subtraction scheme with 
respect to D = 4 — 2e, which in particular entails the familiar coupling renormalization 



2 2e 



g\v)^Z a {g\ V ),e), 



where at the one-loop level 



i-A) 



-i 2 
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+ 0(g 4 
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-N r 



-N F 



;i.i9) 



(1.20) 



1.3 Problems not Addressed 



It may be clarifying to mention two related problems which, while important in their own right, are 
not addressed in this paper. 

(i) It is by no means clear at this time whether (1.7) with (1.8) is the correct action to use in 
continuum QCD. Recent sharpening by Zwanziger [12,13] of Gribov's criticism [14] of the insufficient 
gauge fixing afforded by (1.8) has rendered untenable the convenient prejudice that "Gribov's problem 
may be ignored in perturbative treatments" . It is now quite clear that if one insists on quantizing the 
chromodynamic system with a non-redundant set of degrees of freedom for the gauge field, the ensuing 
restriction of the path integral to a fundamental domain, containing exactly one representative of each 
gauge orbit, affects also the "perturbative" region of small fields. Zwanziger has in fact presented [12] 
a method of approximately "exponentiating" the fundamental-domain restriction, which demonstrates 
the minimum of new effects to be expected from a more complete gauge fixing. His treatment leads 
to the replacement of the action (1.7) by S' = S + "fit, where 7 is a parameter of dimension mass 4 
and of the non-analytic type (1.3), while H is a nonlocal "horizon" functional, 



H = J d D xd D yAZ(x)f Cl 



M-^y)] J dbe A£(y), 

J ca 



(1.21) 
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with M _1 denoting the inverse of the FP operator in the curly bracket of (1.8). (This nonlocal 
term may be replaced with a local one by using path integration over additional auxiliary fields.) The 
drastic low-energy effects produced by this term - in particular, the emergence of the gluon propagator 
of eq. (2.17) below, describing a short-lived elementary excitation of the gluon field - bear several 
intriguing similarities to those emerging from the "extended perturbative" approach discussed here, 
but the precise relationship of the two methods is unknown. What should be emphasized at this time 
is that there is not necessarily a contradiction between them. The Dyson-Schwinger solution for the 
"redundantly" described system of eqs. (1.7/1.8), when given enough nonperturbative freedom, may 
settle down self-consistently in the region singled out by the non-redundant description, and the recent 
insight that the effect of (1.21) can be described in terms of a nonperturbative vacuum spontaneously 
breaking BRS invariance [15] may be seen as pointing in this direction. 

(ii) The present paper studies a purely Euclidean theory. In the few instances where we need 
to refer to its properties in the Minkowskian domain, as in sect. 5 below, we shall proceed, like 
most current investigations of QCD, on the assumption that the usual strategy of Euclidean field 
theory - to define Green's functions and solve for them entirely in the Euclidean, and to continue 
to the Minkowskian only in the final answers - yields physically correct results. For the gauge-fixing 
dependent correlations of the elementary QCD fields, from which any analytic investigation has to 
start, there seems to be at present no full proof of this, and the discussion below will show that in the 
approximating subsequence of primary interest for QCD, the Euclidean solution continued analytically 
will differ from the corresponding direct solution to the Minkowskian equations. Since the core parts 
of the method described here, and in particular the basic quantum effect leading to self-consistency 
of the nonperturbative terms, are general enough to continue to work in e.g. a purely Minkowskian 
theory, we have not given this question any priority, but its existence should be kept in mind when 
proceeding to various applications that the method will invite. 

2 The Extended Perturbative Expansion 
2.1 General Restrictions 

We next describe the form and general properties of the extended iterative sequence, with no attention 
as yet to the question of how this form achieves self-consistency in the DS equations of the theory. In 
this sequence, as the terminology indicates, the organizing principle of the perturbative solution (1.17) 
is not discarded altogether: one still considers a formal power series in the parameter [g(v) / Att] 2 , and 
therefore a weak-coupling solution that is directly applicable (i.e., applicable without infinite resum- 
mations) only if that parameter remains sufficiently small at all scales v to permit semi-convergent 
expansions of the Tjy. It is worth emphasizing that this assumption is entirely compatible with present 
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fragmentary knowledge about the flow of g in QCD: what is truly known of g{is) are a few leading 
orders of its asymptotic expansion at large is, which are of the inverse-logarithmic form of eq. (1.4), 
with k 2 replaced by is 2 . As for low-is behavior, although folklore about "the running coupling blowing 
up at the scale A" has become so pervasive as to be occasionally confused with theory, the emerging 
consensus, if any, from phenomenology [16] and lattice studies [17] seems rather to point in the direc- 
tion anticipated by Gribov in 1987 [18]: as is — > 0, the running coupling does nothing dramatic, but 
"freezes" around an order of magnitude of perhaps 0.2 for the quantity [g(is)/2ir] 2 = a s j-K in the MS 
scheme. In the present context, as we shall see, amplitudes can get large at momenta « A, but via a 
rather different route opened up by the novel feature of the expansion. The new feature is that we 
now have a double, two-index sequence of approximants, in which each term is allowed an additional 
dependence on the g 2 -nonanalytic mass parameter (1.3): 



T N ({k};g 2 (is);is) = hm hm T [ ^\{k}; g 2 (is), is) , (2.1 



*oo p—*oo 
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T^'\{k}-k-is) (2.2) 

P '=i 



(We are using an abbreviated notation suppressing all dependence on parameters not immediately 
relevant to the present argument, such as quark masses and gauge-fixing parameters.) One may 
immediately state two boundary conditions on the nonperturbatively modified rt r ' p l amplitudes. Since 
A 2 , by (1.3), vanishes faster as g 2 — ► than any power of g 2 , it makes physical sense to consider the 
(formal) limit in which A 2 — > but the g 2p remain finite. In this perturbative limit we should ensure 

T [r,p) ( A 2 = ) = r (p)pcrt (p = 0,1,2...)- (2-3) 

On the other hand, since QCD is asymptotically free, and since the logarithmic corrections to asymp- 
totic freedom are known to arise from partial resummation of the g 2 -power series, it is plausible to 
demand that the zeroth-order functions (p = 0) should possess naive asymptotic freedom, i.e. 

({AA:}) r(°) pcrt ({AA:}), A oo , (2.4) 

as the set {k} of external four-momenta are scaled up uniformly. These simple-looking conditions 
will be seen to strongly restrict the nonperturbative extension. We will actually impose a restriction 
slightly stronger than (2.4), namely the requirement that 

the nonperturbative extension of r( pert ) 

y (2.5) 
should continue to be perturbatively renormalizable. 

In an asymptotically free theory, and only there, this requirement is intuitively plausible: the large- 
momentum behavior of vertex functions in the loop integrals of (1.13) is known to be essentially the 
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perturbative one, apart from slowly varying logarithmic modifications that should not change the 
divergence pattern qualitatively. What makes (2.5) a somewhat stronger statement (for vertices with 
N > 3) is the implied condition that behavior no worse than for the perturbative vertex should obtain 
even when only those momenta ofT that run in a loop become large, while the remaining, "external" 
ones are kept constant. Comparison between eqs. (3.8) and (3.9) below, for the three-gluon vertex, 
will pin down this difference more quantitatively in a specific example. 

We already emphasized that in the " nonperturbative direction" of the sequence (2.1), characterized 
by the index r and relating to the dependence on A 2 , the approximation, in contrast to the local one 
provided by a Taylor series around a point, must be global. On the other hand, at the basis of 
any perturbative renormalization process lies the possibility of superficial convergence assessment by 
integer-power counting. The only simple meeting point for these two requirements, and for the direction 
provided by the examples in sect. 1 of nonperturbative quantities, is the idea of representing r[ r, °l, the 
nonperturbatively modified vertex functions of zeroth perturbative order, by rational approximants in 
A 2 of increasing order r. This will indeed be seen to lead to a viable approximation scheme. 

2.2 Gluonic Two-Point Function 

In this section, we give details for the simplest case of a scalar vertex function with only one invariant 
momentum argument: the two-point vertex, or negative-inverse propagator, of transverse gluons. In 
the notation explained in sect. 1, this is 

r ^' 2 » = -i5M- < 2 ' 6 > 

where Dt is the invariant function defined by the tensor decomposition of the full Euclidean gluon 
propagator, 

D^(k) = t^(k)D T (k 2 ) + l^(k)D L {k 2 ) ; (2.7) 
k^k u 

t^(k) = 6^ - = 6^ - l^(k) ; (2.8) 
k z 

D L {k 2 ) = ^. (2.9) 

The propagation characteristics of the gluonic elementary excitation, A^(x)|0 >, will be determined 
by the zeroes and branch points of Tj> in the complex k 2 plane. In the Euclidean domain, its rational 
approximants, which we will characterize by their denominator degrees r, will be of the form 

_ r M],, 2 A2 x (fc 2 )'- +1 +C,4A 2 (fe 2 )-+...+C %T . + i(A 2 )'- +1 , 

1 T [K ,A ) - (fc2)r + ^ )lA 2( fc 2)r-l + ... +f?f . )r (A2)r > 

(2.10) 

r = 0,1,2,3,... ; 
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where the 2r + 1 coefficients Cr,i and rj r j are all real. Note how the boundary condition (2.4), in which 

-rp )pert (k 2 ) = k 2 , (2.H) 



has uniquely fixed both the relative degrees and the leading coefficients of the numerator and denom- 
inator, and how condition (2.3) is then automatically fulfilled. Only in an asymptotically free theory 
does one have such strong a priori restrictions on the form of the approximants. 

Without imposing specific dynamics, the sequence (2.10) still covers a variety of physical situations. 
We do not give a complete classification here, which would include several unphysical or exotic cases, 
but mention just the two subsequences of primary physical interest: 

(1) "Particle" subsequence. Here r is even, so that the propagation function Dp°\k 2 ), by (2.6), has 
odd denominator degree r+ 1, and therefore at least one pole on the real k 2 axis, since the coefficients 
are real. (Note our convention of using the rational-approximation index r also on the propagator, 
although for the latter it gives the numerator, rather than the denominator, degree.) If the real pole 
closest to the origin sits at timelike Minkowskian (i.e. negative Euclidean) k 2 , it represents a stable, 
asymptotically detectable gluon particle. Assuming - as everybody seems to have assumed tacitly 
since the classic papers of Lehmann and Kallen [19] - that the elementary operator field can connect 
the vacuum to at most one single-particle state, one would expect this mass-shell pole position to 

\r 01 

stabilize as r is increased. The remaining r poles and r numerator zeroes of D l T ' would then be 
expected to come again on the real axis, but separated from the particle pole, and for increasing r 
would be expected to settle into the alternating pattern that in the context of rational approximants 
is known [20] to approximate a branch cut - the Lehmann-Kallen dressing cut, arising from virtual 
decays of the particle into multiparticle configurations (Fig. 1A). 
The simplest approximant of this sequence, r = 0, is 

- Tp°\k 2 , A 2 ) = k 2 + Co,iA 2 (Co,i real and > 0) (2.12) 

which when compared to (2.11) represents the general Schwinger mechanism [21]: the spontaneous 
creation of a mass m 2 = Co,iA 2 of the nonperturbative type (1.2) in a massless bare propagator. 

(2) " Quasiparticle" subsequence. Here r is odd, so that Dp ^ has an even number r + 1 of poles, 
and at least one real zero. Ignoring again the exotic possibility of two or more stable-gluon poles 
at different real masses, one would expect the two poles closest to the origin to come as a complex- 
conjugate pair at, say, k 2 = — <7 rj ±A 2 , with <r r>+ = a*_ being a dimensionless complex number. As 
discussed in [22] and in sect. 5 below, the conceptual problems apparently caused by this complex- 
poles structure are not insurmountable, provided the solution is used consistently. The pair would 
represent an intrinsically short-lived elementary excitation of the gauge-vector field, with lifetime of 
the order of 1/A. Note that the real zero closest to the origin - corresponding to a pole of the function 
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(2.6), i.e., a singular gluonic self-energy - is every bit as essential in this context as the complex pole 
pair, as it will be seen to make the connected Green's functions nonsingular in the invariant masses 
of external vector lines. The remaining even number r — 1 both of poles and of zeroes may then come 
in complex-conjugate pairs, farther from the origin than the leading " quasiparticle" pair, with poles 
and zeroes interspersed so as to approximate two conjugate branch lines (Fig. IB). These again would 
represent dressing - the short-lived excitation coupling dynamically to multiple copies of itself. The 
absence of real-axis branchpoints would signal that it has no channels for decay into stable fragments 
that would jointly carry the open quantum numbers of a gluon. This structure would fit most closely 
the empirical situation for the short-lived gluon presumed to be present at the origin of a gluon-jet 
event. 

The simplest approximant of this sequence, r = 1, is 

(A; 2 ) 2 + Ci,iA 2 fe 2 + Ci, 2 A 4 



-rg >0I (fc 2 ,A 2 ) 



with real coefficients given by 

= Ci,i - mx> 

and satisfying 



k 2 + r/i iA 2 



= k 2 + «, iA 2 + 



Ul,3 



A 4 



k 2 + m 2A 2 



(2.13) 



ui,2 = mX' u i,3 = Ci,2 - »7i,i«i,i; 



Ul,3 > 



;(«l,l-«l,2) 



The corresponding nonperturbatively modified propagator of zeroth perturbative order, 



D^ 0] (k 2 ,A 2 ) 



k 2 + Ui 2 A 2 



(2.14) 



(2.15) 



(fc 2 + ( j 1 , 1 A 2 )(£: 2 + cT 1 ,3A 2 ) ' 

exhibits the minimum of features mentioned above: a real zero at k 2 = — «i 2A 2 , and a "quasiparticle" 
pair of complex-conjugate poles at k 2 = — cr^A , where 



= 0-1,1 = j ("1,1 + "1,2) + iJ "1,3 



1, 

2K1 



1 2 



"1,2, 



^1,3 = a l - 



(2.16) 



At r = 1, there are no further zeroes and poles as yet to simulate dressing cuts. 

Eqs. (2.15) and (2.13) are of a form suggested in [22] (for u\p = 0) and used in [7] as an 
element of an approximate, nonperturbative DS solution, and that form can now be identified as the 
lowest member of (the odd-r subsequence of) a systematic sequence (2.10) of approximants to the 
nonperturbatively modified function (2.6). A special case, 

k 2 



D 



(GZ) 



(k 2 ) = 



fc 4 + 7 ' 



(2.17) 
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with ui : i = ni 5 2 = and ui^A 4 = 7, had actually been arrived at much earlier by Gribov [14], and was 
later derived independently by Zwanziger [12], via the entirely different route of fundamental-domain 
restrictions. Because of this different origin, the Gribov-Zwanziger 7 term is present already at tree 
level, i.e. in the analog of what we called r(°) pert , whereas the mechanism to be used in sect. 4 to 
stabilize nonzero coefficients in (2.13) will be seen to operate from DS loops. 

It is clearly desirable to try out all of the above types of approximants in the DS equations, to 
determine those that can achieve dynamical self-consistency and, if necessary, to further distinguish 
between the latter by a stability analysis. We do not here embark on such a comprehensive study, 
but describe a few steps toward the much more restricted program of trying out the " quasiparticle" 
subsequence - with odd r and only complex-conjugate propagator singularities. In this subsequence, 
the nonperturbatively modified two- vector vertex in zeroth perturbative order will be of the form 

u r , 3 A 4 



-Tp 0] (k 2 ,A 2 ) 



k 2 + u r 1A 2 + 



k 2 + n ri2 A 2 



(r-l)/2 

+ E 

s=l 



U r ,4s+l 



A 4 



+ 



A 4 



k 2 + n rj4s A 2 k 2 + n ri4s+2 A 2 



(2.18) 



r = 1,3,5... 
generalizing eq. (2.13), with u r ,i, Ur,2> ^r,3 real, with 

U r ,4s+2 = U*,4s' ^r,4s+3 = K,4s+1 



8=1... 



r — 1 



and such that all poles of the corresponding propagation function, 

k 2 + u r 2 A 2 



Dp 0] (k 2 ) 



(r ^ 2 (k 2 + u rM k 2 ){k 2 + u rM+2 K 2 ) 



(2.19) 



(2.20) 



(2.21) 



(k 2 + a rA A 2 )(k 2 + a r , 3 A 2 ) ^ (k 2 + a r , 4s+1 A 2 )(k 2 + a rAs+3 A 2 ) ' 

form complex-conjugate pairs, with a rj i,a r ^ denoting the leading "quasiparticle" pair closest to the 
origin of the k 2 plane, i.e., 

r — 1 , 



°V,3 — f*,i; Cr,4s+3 — Cr*,4s+l'i l°V-,4s+l| > l°V,l| (s — 1 . . . 



(2.22) 



The naive-asymptotic-freedom condition built into (2.10), Tp°^ — ► — k 2 for k 2 » A 2 , guarantees that 
in the pole decomposition, 



Pr,2 



k 2 + a rl A 2 k 2 + a r3 A 2 



(r-l)/2 

+ E 



s=l 



Pr,4s 



+ 



Pr,4s+2 



k 2 + fT ri4s+ iA 2 k 2 + OY i4s+ 3A 2 



the sum of the dimensionless residues is unity: 

(r-l)/2 



Pr,0 + Pr,2+ (Pr,4s + />r,4s+2) = 1 (r = 1, 3, 5 . . .) . 



(2.23) 



(2.24) 



s=l 
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This, incidentally, is one of the features that distinguish the functions (2.23) from propagators plagued 
by so-called ghost poles, which - usually as a result of inadequate approximations - vexed field theorists 
in the 1950's [23], and for which the sum-of-residues was negative or vanishing. 



2.3 Relation with the OPE 

To exhibit the connections of the sequential approximation (2.10) with an established area of QCD 
methodology, we briefly look at the operator-product expansion (OPE). The OPE for the elementary 
QCD fields has only recently begun to be established correctly (for recent results and literature, see 
refs. [24] through [27]) and is not commonly discussed in the present terms, yet it already represents a 
step towards describing the "truly nonperturbative" A 2 dependence. Again we consider the simplest 
case of the one- variable vertex function (2.6), and restrict ourselves to a theory with at most massless 
quarks, so as not to have to worry about the role of other invariant mass scales besides A. By writing 
an OPE for the Euclidean two-point function of the gauge field, transforming to momentum space, 
contracting with the #*" of eq. (2.8) to project out the transverse portion, and forming the inverse 
function (2.6), one arrives at an expansion of the general form 



r T (k 2 ;g 2 (v);v) 




k 2 



+A 2 L 10 + E 



I p=l 



(2.25) 



+ ... 



(A 2 ) 



'W) 



+ ... 

where the coefficient functions (modified Wilson coefficients) are p — th degree polynomials of loga- 
rithms: 



L _ JO) , (1) j f ^_ 



+ 



+ C 



(P) 
n,p 



In 



k 



(2.26) 



This may be viewed as a form of the series (2.1/2.2) where the r[ r ' p ) have in turn been expanded in 
power series in A 2 . More precisely, apart from the typical perturbative logarithms, it is an expansion 
(presumably semi-convergent) in powers of (A 2 /k 2 ) for k 2 » A 2 , the deep-Euclidean limit. The 
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first line of (2.25) is identical with the perturbation series (1.17), so the boundary condition (2.3) 
is satisfied. The presence of the other terms, containing powers of the g 2 -nonanalytic scale (1.2), 
shows clearly that the perturbative series alone would be an incomplete solution even when summed 
to all orders. At the core of these additional terms are the quantities L n o(A 2 )™,n > 1, which the 
OPE derivation identifies as linear combinations of vacuum expectation values of composite operators 
("vacuum condensates") of increasing mass dimension, 2n. 

As it stands, the OPE (2.25) does not satisfy our needs, for two reasons. First, the OPE by 
itself does not determine the quantities L h q - this requires a truly dynamical principle, such as the DS 
equations. Second, even if the L„o were determined dynamically up to some n, no finite order of (2.25) 
would be a satisfactory continuation of the vertex function into the region of primary physical interest 
- the region k 2 < A 2 of typical hadronic masses. In this region, given our assumption that [g(i/)/<iir] 2 
never becomes large, the important task clearly is to obtain a continuation-through-resummation of 
the "vertical" sums in (2.25) - in particular, the p = vertical sum, 

-Tp 0] (k 2 ,A 2 ) = k 2 

which is free of perturbative logarithms. (Once this is done, the remaining (p > 1) "vertical" sum- 
mations can in principle be generated iteratively from r^' ! through the DS equations, just as in the 
perturbative case.) From this standpoint, the rational approximations (2.10) may now be viewed as a 
systematic sequence of continuations-through-resummation of the OPE subseries (2.27), which locate 
its \ow-k 2 zeroes and singularities with increasing accuracy. 

We stress that the OPE (2.25) has been used here only as a point of comparison. The "extended- 
iterative" scheme differs from it in more than the technical aspect of parametrization. (In the OPE, the 
fundamental parameters are an infinite set of dimensionful vacuum condensates, whereas in (2.10) that 
role is played by the dimensionless vertex coefficients such as (Cr,i, Vr,j), and condensates are secondary 
quantities calculable in principle in terms of the latter - see the second of refs. [7] for [r,p] = [1,0] 
examples). Its properties can be qualitatively different from any finite order of (2.25) because the two 
are separated by a nontrivial step of analytic continuation. In particular, the "horizontal" sums in 
(2.25) are always ordinary QCD perturbation series based on the Feynman rules, with all the attendant 
problems arising from the empirically wrong zeroth-order spectrum (including free, massless, physical 
gluons, and their associated infrared singularities) of the latter. By contrast, expansion (2.1/2.2) is 
based on the use of (inter alia) modified propagators such as (2.21), which in particular produce no 
infrared singularities at all (the unphysical longitudinal gluons still do, but their effects are always 
preceded by £ factors that identify them as gauge-fixing artefacts). This opens up the possibility 
that Borel transforms with respect to g 2 of eq. (2.2) may have, apart from gauge-fixing artefacts, no 
genuine infrared renormalons, the effects of the latter having been absorbed in a redefinition of the 



A 2 /A 2> 



+ 



(2.27) 
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zeroth perturbative order. 

3 Three- Vector- Vertex Approximants 

The hierarchical structure of the DS equations (1.13) implies that in principle all proper vertices 
Tjy should be treated simultaneously by mutually consistent, nonperturbative approximants. The 
method described here will however be found to have the simplifying feature that the essential self- 
consistency problem - for the vertices T^' " 1 - is rigorously restricted to the small finite set (1.14) of 
superficially divergent vertices. In this section we therefore compile formulas for rational functions 
~^ sdiv concentrating on the example of the proper three- vector vertex I^y, which serves to illustrate 
all the essential features. The largely analogous material for the remaining T S( n v vertices will be 
relegated to appendix A. The only new aspect this discussion will turn up will be the " factorizing- 
denominator" rule discussed in connection with eq. (3.6) below. Otherwise, the main complications 
will be the multi-variable nature, and the notoriously unwieldy color-and-Lorentz-tensor structure, of 
the N > 3 QCD vertices. One measure we will adopt to control this purely kinematical complexity 
is to restrict ourselves to vertices with at most transverse external gluon lines, wich are sufficient for 
performing calculations in the Landau (£q = 0) gauge. As shown most clearly by the example of eq. 
(2.9), it is in these that nonperturbative effects are expected to develop most freely. Amplitudes with 
at least one longitudinal gluon are strongly restricted by the Slavnov- Taylor (ST) identities, whose 
nonperturbative saturation is not a subject of this paper. 
The three-gluon vertex has color structure, 

in terms of f-type (antisymmetric) and d-type (symmetric) structure constants. Each color component 
in turn is of the general Lorentz structure discussed by Ball and Chiu [28]: a linear combination of 
14 independent third-rank tensors, with 6 different invariant functions having appropriate symmetry 
or antisymmetry properties in the Lorentz-scalar variables p\,p\-,p\- We consider only the totally 
transverse portions, 

r 3T"l) = t^(pi)t x '\p 2 )t^(p 3 )T^ {c) (c = f or d) , (3.2) 

which have contributions from only 4 Lorentz tensors combined with two different invariant functions 
Fq,Fi for each color component: 

r 3T A ('c)'(^^2,P3) = ^^p 1 )t X '\p 2 )t^(p 3 )x 

*{5 Xu (P2-PzYF (c] ,(pIpIp\) 
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+ ^( P3 - P i) X F {c)0 (pIpIp 2 2) (3-3) 
+^\pi-P2) u F {c)0 (plplpl) 
HP2-PzY{pz-Pi) X {pi-P2) u F (c)1 {pIpIpI)) ; (c = / or d) . 

Remember pi + P2 + Pz = 0. Here the dimensionless functions Fr c \ are symmetric (for c = /) or 
antisymmetric (for c = d) in their first two arguments, while the functions ^(c)i> °f mass dimension 
—2, are totally symmetric (c = /) or antisymmetric (c = d) in all three arguments. The perturbative 
zeroth-order limits are, 

F (c)k eit = M*o (c = /ord, fc = 0orl). (3.4) 

In setting up sequences of rational approximants for these four invariant functions, capable of dynami- 
cal consistency with the gluon-propagator sequence (2.18), one encounters a new aspect: approximants 
with the most general denominator polynomials in the three variables p\,p\iP\ are not useful. The 
zeroes of such a general denominator in any one variable pf are complicated non-rational functions of 
the two other variables, and this will stand in the way of DS self-consistency when a vertex transfers 
its structure in an external momentum to the next lower vertex through the hierarchical coupling. At 
the expense of slower convergence of the approximating sequences, we must restrict ourselves to the 
narrower but still sufficiently general class of factorizing- denominator rational approximants (FDRA), 
i.e. those in which the denominator factorizes with constant zeroes in all three variables. To see 
that such a more special approximation is always possible in principle, consider e.g. the function 
F(f)o(PiiPhp1) as the kernel of a symmetric integral operator parametrically dependent on p§, and 
write an eigenfunction expansion 

F U )0{plphpl) =Y,9n{pl)fn{pl)fn{pl) , (3.5) 
n 

with eigenvalues g n . By using single- variable rational approximation for the f n and g n functions, 
letting n range over finite but increasing numbers of eigenvalues, suitably discretizing the integral 
in case there is a continuous spectrum, and putting everything over a common denominator, one 
generates a sequence of approximants in which denominators have the desired, fully factorized form. 
Similar considerations apply to the other invariant functions. 

By extension of our above definition, the degree r of rational approximation for the three-point 
vertex in zeroth perturbative order will be the number of different denominator zeroes in any one 
variable for the entire tensorial vertex. That is, to order p = in eq. (2.1) and degree r, all invariant 
functions will be of the form 

„[r,0], 2 2 2n N 3V(c)k(-P^P^Pl) . . 

F[ c)k J (Pi,P 2 ,P 3 ) = TT— ; i r r ; 7TTT— 7T ( 3 -6) 



n (pi + < 2s a 2 )J [n (p\ + < 2s a 2 )J [n (*§ + < 2s a 2 ) 
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(c = / or d; k = or 1; r = 1, 3, 5, . . .) 

with the same fully factorized denominator, where the appropriate mass dimensions and Bose - symme- 
try restrictions, as well as the boundary conditions (2.3/2.4), are built into the numerator polynomials 

<W (A Arf) = M» (*«)'+ E 4&» (*>?)"" (ri)"" (*>§)"" ( A f 

mi,m2,rri3>0 

(3.7) 

Condition (2.4) allows nonzero coefficients C mim2m . A only for 

m\ + m<i + m3 < 3r — — 1 . (3.8) 

Additional restrictions follow from the postulate (2.5) of preservation of perturbative renormalizability. 
As a minimum, it requires that no part of Y^y should lead to ultraviolet divergences stronger than the 
corresponding perturbative ones. Since in a 1PI diagram one of the three legs of T^y may be external, 
we conclude that when any two of its momenta are running in a loop, the vertex should behave no 
worse than q l at large loop momenta q. For the invariant functions (3.6) this requires 

mi + rri2 < 2r — k, mi + m% < 2r — k, + m\ < 2r — k . (3.9) 

Moreover, the coefficients should exhibit the symmetries with respect to mi, mi,mj, necessary to fulfill 
the conditions of partial or total (anti-)symmetry listed before. 

Comparison between eqs. (3.9) and (3.8) illustrates well the statements made after eq. (2.5): of 
the conditions required for a general rational function to reduce to one that preserves perturbative 
power counting, the vast majority are already enforced by the asymptotic-freedom condition (3.8), 
with (3.9) representing only a mild additional restriction. (In fact, on the r = 1 level written out 
below, (3.9) will already be fully implied by (3.8).) For the discussion of sect. 4, the essential and 
simple consequence of the factorizing-denominator structure is that with respect to any one argument 
pi, the full p = vertex has a partial-fraction decomposition 

^<™(f>l,i> 2 ,^> 3 )f'' , = B^feft*) 



s=l 



k = 1, 2, 3 and cyclic) , 

where invariant functions for the residue tensors B s ^ r with s > 1 depend only on pf and p|, whereas 
the regular part i?o,r is in addition an r-th order polynomial in p\. (Note that with respect to any 
single pi, the conditions (3.8/3.9) do not rule out terms that grow polynomially. Without dwelling on 
this point, we remark that most of these terms will however turn out to vanish on dynamical grounds.) 
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The question may be raised as to whether the choice of variables p\ , p\ , p 3 in constructing the FDRA 
of eq. (3.6) is unique. Could other sets of Lorentz invariants be used, such as {p± • p 2 ,p 2 • P3,P3 ■ Pi}? 
In eq. (3.6), we have, for brevity, anticipated the fact that dynamical consistency between the Y 2 v 
and Fsv vertices will indeed require that FDRA's for r 3 y be constructed in the pi variables. 

To illustrate these structures, we write expressions (3.6), now fully decomposed into partial frac- 
tions in analogy with eq. (2.18), for the simplest case, r = 1. In this case, all rational structure can 
be expressed in terms of the single-variable pole factors 

A 2 



IT 



(i = 1,2,3) 



Pi + <2 A2 

The f-type invariant-function approximants then read, 

*(/)o (pIpIpD = 1 + (ni + n 2 ) + x li3 n 3 + (x 1>2 + x' lj2 ±) m 2 



(3.11) 



+ 



(x 1A + x' 1A ^) ni + (x 1A + x[ A ^j n 2 



n 3 



(3.12) 



+si )5 (M 2 n 3 ) , 

^(/fi (pIpIpI) = ^2 [ x i.6 (nin 2 + n 2 n 3 + n 3 no + Xl , 7 (Mails)] 

whereas the d-type approximants are, 

^K 1 {pIpIpI) = xi,8 - n 2 ) + ^!, 9 + x' li9 i-) ni - ^ li9 + si )9 i-) n 2 



n 3 



[1.0] 



0. 



(3.13) 



(3.14) 



(3.15) 



(d)i \PnP21P3j 

The latter result reflects the fact that for k = 1 and r = 1, no fully antisymmetric numerator 
polynomial in three variables obeying restriction (3.8) can be constructed. 
To order (r = l,p = 0), the set of twelve dimensionless coefficients, 



{xi,i,xi, 2 , . ..x[ 9 } , 



(3.16) 



which must be real numbers of order unity, are the fundamental vertex parameters to be determined 
from the DS equations (1.13). Note that all invariant functions at this level feature one single pole 
position, — n' 12 A 2 . While e.g. in the k = functions the partial (anti-) symmetry would seem to 
permit the use of two different pole positions, such an approximant, under our above definition, would 
already be of degree r = 2. This strict definition of r may at first appear to be overly rigid, but it will 
turn out to be the relevant one for the self-consistency problem. 
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4 Dyson- Schwinger self-consistency 




When introduced into the DS equations, the nonperturbative sequence (2.1), in contrast to the purely 
perturbative one, encounters a nontrivial self-consistency problem - the self-reproduction of the p = 
functions r[ r '°l at a given level r. For the first iteration, and in the notation of eq. (1.13), this problem 
can be stated as 

r [r,0]] 1 = [pM) _ r (0)pertj + Q ^2 {v) <r + ^ (4J) 
J R,u 

Here, e(r + 1) stands for the approximation error in the r direction which, in contrast to that of the 
semiconvergent p sequence, cannot be characterized by a power of some expansion parameter. 

It is useful to reflect at this point on what the precise meaning can be of "solving a Dyson-Schwinger 
system by a sequence T^ r ' p ^: 

(1) Upon iteration of the system around the starting point r[ r '°l, the sequence should reproduce 
itself at each iteration step up to corrections of the next higher perturbative order p, which are simul- 
taneously generated in the process. 

(2) Self-reproduction of rj^' ', which features a number n(r) of dimensionless nonperturbative 
parameters, can only mean that this approximant is made to coincide with the r.h.s. of its own DS 
equation at n(r) points in the space of scalar momentum variables k 2 of IV, or more generally (since 
global approximants allow a variety of matching prescriptions) with respect to n(r) "comparison data" 
(function values, derivatives, pole positions, integrated quantities ...). At all other points, or in all 
other data, a matching error will necessarily remain that can be improved only by going to the next 
higher r. 



4.1 Transverse-Gluon Self-Energy 

It is characteristic of the scheme discussed here that the hardest problem for direct numerical Dyson- 
Schwinger solutions - the self-reproduction of momentum structure - has a relatively straightforward 
and explicit answer: the (N' > iV)-point vertices entering the dressing functional <& will "hand down" 
their rational structure in external momenta to the N-point ones. This pattern may at first sight 
appear simplistic, but in fact turns out to be an efficient way of exploiting the most difficult and 
peculiar structural feature of a Dyson-Schwinger system - the hierarchical coupling. It has, however, 
the immediate consequence that the set of nonperturbatively modified vertices can only be treated 
as a whole: no self-consistency, however approximate, will be possible in this framework if one seeks 
solutions to particular DS equations while treating the higher vertices appearing there by more or less 
unrelated assumptions. 
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Less straightforward is self-reproduction of the coupling structure. This problem is already clearly 
visible in (4.1): the DS loop integrals constituting <I>jv must produce the bracket on the r.h.s., with no 
g 2 prefactor, in spite of the fact that they always come with at least one g 2 , prefactor. The importance 
of securing this feature can hardly be overstated: since perturbative corrections of finite order (p > 1) 
cannot alter the qualitative spectral properties of a solution, the essential nonperturbative features, 
and in particular qualitative changes expected in the spectrum of elementary excitations - i.e., in the 
two-point functions of the basic fields - must establish themselves already on the p = level. 

The mathematical mechanism leading to the " eating of g 2 prefactors" , a low-order version of which 
was used already in [22] and discussed in detail in [7], will again be demonstrated for the example of 
the one- variable function (2.6). Its DS equation, 



-k 2 + 



-V E 

An) z - 



4? [r 2 ,r 3 , r VVVT ] , 



(4.2) 



r T (fc 2 ) 

M=A...F 

is stated diagrammatically in fig. 2, showing the one-DS-loop (M = A, B, C, D) and two-DS-loop 
terms (M = E, F) of the $tt functional. (The term "DS loop" will be used to refer to the specific 
structure of DS integrals which are neither bare nor fully dressed, containing always one bare vertex 
times a number of dressed functions). Start at the one- loop level, where the terms M = E, F do not 
yet contribute. For eq. (4.1), we have to evaluate 



<I> 



(l=i) 

TT 



p[r,0] 



<l> 



(A) 
TT 



D [r,0] r Ml 
1J ' VTV 



(B) 
TT 



+4? 



+ E 4? 



f)[r,0] Arfi] 
' GTG 



,[r,0] -.[r.O] 



(4.3) 



FTF 



Now use the FDRA structure of the [r, 0] three-point vertices, as expressed in the partial-fraction 
decompositions (3.10), (A. 24), and (A. 44). Expression (4.3) then decomposes into terms regular and 
singular with respect to the external variable k 2 : 

A 2 \ 



(i=i) 

TT 



r Mi =ij[\k 2 ) + Y / 4 r) (k 2 



s=l 



k 2 + 



r,2s 



A 2 



(4.4) 



The quadratically divergent loop integrals /i^ are of the same form as (4.3), but with the replacements 



7,0] 
VTV 

r,0]_ 
GTG 

r,0]_ . 
FTF 



B (r) 


in 


^TT ' 


B% 

(f)s 


in 


^TT ' 


Mr) 
C (F)s 


in 


TT ■ 



(4.5) 

(4.6) 
(4.7) 



Indices (/) on the r.h. sides of (4.5/4.6) recall the fact that only the f a b c parts of those vertices 

(c) 

contribute here. By definition, the "tadpole" contribution (which does not vanish as in the 

perturbative case), is entirely included in Iq T \ as it has no /c 2 -singular terms, and in fact is a constant. 
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In eq. (4.4), since we have treated the 3-point vertices consistently, i.e. by FDRA's of the same 

(r) 

level r, the r.h.s. already displays the needed number of poles, and inspection of the Is integrals 
shows that they produce logarithmic branch points (generally complex since we are using the complex- 
pole propagators of the odd-r sequence) but no poles in k 2 of their own. It is then natural to choose 
the r positions and r residues of those poles as comparison data - if the poles were not matched, the 
matching error would be locally infinite. That is, one requires 

u r,2s = u' r . :2s {s = 1 . . . r) ; (4.8) 

- u r , 2s+1 A 2 = 2 IP (-u r , 2s A 2 ) (s = 1 . . . r) . (4.9) 

Eq. (4.8) expresses the "handing down" of momentum structure, rationally approximated, from the 
3-point vertices to the 2-point one under consideration. It immediately implies, by eq. (2.21), that 
the propagator will have zeroes at the positions of the gluon-variable poles of the 3-point vertices. 

It would seem that as a (2r + l)-th comparison datum to fix the one still undetermined coefficient 
u rt i of (2.18), we could use the value of the smooth-remainder function, 



tfV) = Io\k 2 ) + A 2 ± ^^)-I ( l\-u r , 2s K 2 

s=l 

at some arbitrary point, k 2 = —uA 2 , of the k 2 plane 



jfV)=W + A*S * ">-'• (4.10) 

s=1 K ~r «r,2s ii - 



u r ,,K 2 = {^f 4\-uK 2 ) . (4.11) 



We will soon see that the arbitrariness introduced at this point is only apparent. In eqs. (4.9) and 
(4.11), the problem of how to "eat" the <? 2 prefactor to produce a p = quantity is still with us. 

4.2 The self-consistency mechanism 

Evaluating now the I s and Jo integrals in dimensional regularization - note that all these integrals, 
because of the rational structure of integrands, can be evaluated by standard methods -, one finds 

(s = l...r), 
(4.12) 

J^\-uA 2 ) = (^-J^j ^a (u,v,w,x,y,z) • ^ + [terms finite as e -> 0] + 0(e) j- A 2 , (4.13) 

where the notation recalls that the A scale, too, should be continued to D = 4 — 2e through the 
replacement f3(g') — > P(g') — eg' in its definition (1.2). The sets u. . . z are the dimensionless-coefficient 
sets of the [r, 0] propagators and 3-point vertices, as defined in sect. 3 and appendix A, and include, 
for the time being, the quantity u of eq. (4.11). 



/j r ^(— n rj 2 S A 2 ) = ( — j \ |a s (u, v, w, x, y, z) ■ - + [terms finite as e — > 0] + 0(e) | A 2 
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To analyze the product (go/47r) 2 (A 2 /i/ 2 )~ e (l/e) now appearing in the self-consistency equations, 
we next impose the condition, inherent in (2.5), that perturbative coupling-constant renormalization, 
eq. (1.19), should apply. Then this product becomes, 
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(4.14) 



To get a feeling for its structure, we first use a slightly sloppy argument [7,6] which, however, happens 
to convey the essential points. Start from 
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-f = l-eln ^ +0(e 2 ), 



which gives 



^] +- + 0(e). 



(4.15) 



(4.16) 



The logarithmic term survives the removal-of-regulator limit, e — ► 0, in the same, familiar way as do 
the perturbative logarithms of, say, eq. (2.25), but it differs from them in one important respect - it 
leads them by one order in the perturbative g 2p classification. From eq. (1.3), at the one-loop level 
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This is precisely the 1/g 2 factor needed to "eat" the overall g 2 factor in (4.14); it has survived due to 
its association with the l/e divergence. We now have 
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(4.18) 



But Z a , to the same order, is given by (1.20), so that to the order calculated, the divergences actually 
cancel in the product (4.14): 
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1 + 0(g 2 In g 2 )+0(g 4 )+0(e) 



(4.19) 



Although the terms denoted 0(g 4 ) do include ( as do those in (1.20) ) terms of the form g 4 /e, they 
are not to be kept in a one-loop calculation. We therefore find that not only do the nonperturbative 
parts of (2.18) - the r pole terms and the mass-type constant term - establish themselves in order 
p = in spite of the p = 1 prefactor, but they do so without divergences to the order calculated, so 
that in particular no nonlocal counterterms are needed for the pole terms. This crucial prerequisite of 
perturbative renormalization is therefore automatically preserved. 

The sloppiness of the above derivation based on 1-loop results lies hidden in the terms summarily 
denoted "0(e 2 )" in (4.15) and "0(g 4 )" in (1.20). When examining these more closely, one realizes 



22 



that (A 2 /V 2 )~ e actually has an infinite subseries of terms of type (e/g 2 ) m ,m > 1, whereas Z a has a 
subseries of terms of type (g 2 /e) n , n > 0, and the desired terms of type e (<7 2 ) in (4.14) can therefore be 
produced in an infinite number of ways. To account for all of these, one could obtain both subseries 
exactly by resummations of the LL type, but the result, while perfectly sufficient for the one-loop 
calculation, would not display the basic simplicity of the situation: the quantity II of (4.14) is actually 
independent of g 2 . This follows immediately from the exact integral representations of the two main 
factors. From the extension of definition (1.2) to e / 0, one finds by simple manipulations that 
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where k = [g(^)/4-7r] 2 , and 
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The 1/k factor of (4.20) is the exact counterpart of the one of (4.17) and is what is needed for the 
"eating" process. The lower integration limit m, fixed conventionally as part of the definition of the 
renormalization scheme, will in general depend on e, 



ki =/ei(e) = Ki(0)[l + O(e)], 



(4.22) 



with Ki(0) having the meaning of a renormalized coupling at the scale A. On the other hand, differ- 
entiating (1.19) with respect to v leads to the relation 

d . _ , , , , 1 



K -lnZ Q (.H, e) = - TT|7R , 



(4.23) 



whose reintegration with respect to n under the initial condition Z a (0,e) = 1 gives 
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Z a (n(is),e) = exp 



dX 



A + e/(A) 



(4.24) 



a representation known to t'Hooft [29] in 1973. Now the dependence on n{y) cancels exactly in the 
product of (4.20) and (4.24); the result 

Ki(e) 
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dX 



X + e/(A) 



(4.25) 



has none of the apparent g 2 corrections of (4.19). That 11(e) is, moreover, finite as e —> - the exact 
counterpart of the divergence cancellation observed on the way from (4.18) to (4.19) - follows by 
writing 

IT7M = ITWA) +eft(A ' £) ' (426) 
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where f?i(0, e) = 0; then 
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(4.27) 



(4.28) 



Further subtraction of £>i at A = shows that the integral develops also a term In e, so finally, 

11(e) = -j-[l + 0(e,e In e)] . 
Po 

This is a compact formulation of the self-consistency mechanism. It clearly shows, through the \ 
factor in (4.14), that the mechanism is tied to the divergence structure of the theory, and represents a 
dynamical exploitation of that structure. In this it is reminiscent, of course, of quantum anomalies, 
except that here the process seems to have no connection with the breaking of a classical symmetry. 
The result is that by imposing the algebraic self-consistency conditions 
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■a s (u,v,w,x,y,z) (s = 0...r) 



(4.29) 



in addition to (4.8), the nonperturbative terms of rjp ' reproduce themselves not only without di- 
vergences but also without finite 0{g 2 ) corrections. Note that (4.28) contains only the scheme- 
independent, leading /3-function coefficient, and that it requires evaluation only of the divergent parts, 
a s , of (4.12/4.13) - the finite parts of those quantities only give terms 0(e) relative to (4.29). 

To illustrate these relations, we write conditions (4.29) for the simplest nontrivial, odd-r case 
of r = 1, in Landau gauge, and for massless quarks, where eqs. (A. 28 - 30) apply with mp = 0. 
They are obtained by evaluating the one- loop quantity (4.3) with the r = 1 vertices and propagators 
of section 3 and appendix A, using entirely standard techniques of dimensional continuation, Feyn- 
man parametrization, and symmetric integration, and by extracting the divergent parts as in eqs. 
(4.12/4.13). One finds 
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+5xi,i - gXi >3 + y xi, 2 + 6x M - jXl, 6 



+ ^1,2(21/1,2 + 2/1,5 + yi,9) + jUl^l^} 
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(4.30) 
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(4.31) 



plus, of course, the u' 12 = u i,2 of eq. (4.8). (An obvious simplified notation has been employed for 
the fermion- vertex coefficients z of (A. 45/46).) Together with analogous conditions for the ghost and 
quark self-energies, these will form an algebraic system (linear at this level) determining u\ . . . 1/3, and 
the other self-energy constants, in terms of the nonperturbative coefficients x, y, z of the three-point 
vertices. 

After this extraction of p = nonperturbative parts, the remainder of the DS dressing term, 
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still has a logarithmic divergence. In this remainder, the quantity (4.15/20) no more occurs ( the 
noninteger power that does occur contains the scale k 2 ), and treatment of its coupling structure, 
under boundary conditions (2.3) and (2.5), must therefore revert to the "normal" , perturbative pattern. 
Isolating the divergent piece with R = MS gives, 
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(4.33) 



which apart from the u term is precisely the ( Landau- gauge ) perturbative one-loop divergence. This is 
as expected, since by the very construction of our vertices, (4.32) must contain the correct perturbative 
limit. Appealing now once more to our basic condition (2.5) to demand that this divergence should 
be removable by the perturbative counterterm, one finds that the choice 



u = 



(4.34) 



is dictated uniquely, eliminating the apparent arbitrariness in the matching condition (4.11). This 
observation is particularly interesting in the (non-QCD) context of the even-r "particle" sequence 
with a real-axis propagator pole, since it shows that a nonzero mass term, u r ,iA 2 , can in principle be 
generated while using only the massless perturbative counterterm. 
The finite term of (4.32) in R = MS, 
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where 



i4\k>)- J ° ){k2) - J ° ){0) 



(4.36) 



represents the renormalized, one-loop, radiative correction, having errors both of order g 4 in the 
perturbative and of type e(r + 1) in the nonperturbative direction. 
A number of comments are in order here: 

(i) The divergence content of the one-loop dressing functional is now exhausted. With [r, 0] input, 

it has just been sufficient for producing the divergence (4.33) in the perturbative correction, 
curable by perturbative means, and for triggering the "eating" mechanism of eq. (4.28) in the 
r + 1 nonperturbative terms, as needed for self-reproduction of the input. It does not, therefore, 
generate the terms of the next higher order (r + 1) automatically, as it does in the perturbative 
direction. Even less can it define a unique splitting, valid uniformly at all k 2 , into order— g° and 
order— g 2 terms; this can only be approached asymptotically in the r » 1 limit. At finite r, 
output-input matching can be achieved only in a finite number of comparison data. 

(ii) While perturbative renormalization has now gone through completely at one loop, one notes that it 

has done so with a peculiar new twist: for the vertex function as a whole, the renormalization is no 
more multiplicative, as in pure perturbation theory, since the nonperturbative terms, by (4.28), 
have established themselves in a finite manner. At best, one may write the renormalization 
process as 



with Z3 the perturbative gluon-field renormalization constant, but not as one overall rescaling. 
(This property must hold already in the OPE if one wants to avoid nonlocal counterterms.) 

(iii) The observation that the self-consistency mechanism is tied to the divergences of DS loops 
immediately implies that nonperturbative terms of the present type cannot establish themselves 
in superficially convergent vertices at one loop. The argument can be extended, however [7]: an 
n-loop contribution to a superficially convergent vertex has prefactor (<7o) n but can have at most 
n — 1 divergent subintegrations, which through eq. (4.28) can "eat" at most n — 1 powers of g 2 ,, so 
that the result is at best an ordei-g 2 correction but never of order g°. This is noteworthy because 
it shows that, at least for the limited purpose of determining zeroth-order nonperturbative parts, 
the infinitely coupled nature of the DS system can be beaten without introducing "decoupling" 
approximations: the self-consistency problem of T^' ! is strictly confined to the set of seven DS 
equations for the superficially divergent vertices. 

(iv) The fact that the one-loop evaluation of $ has produced both terms of order g° and g 2 demon- 
strates that there is a decoupling of the loop and perturbative orders, I and p - again a feature 




(4.37) 
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presumably present in any truly nonperturbative solution. A discussion of this aspect has been 
given in the first of refs. [7] and need not be repeated here. Strictly speaking it implies that 
to fully characterize an approximant T^- r ' p \ one would need still another index /, of a technical 
nature, indicating from what loop order I > p the r.h. sides of the 0-th order self-consistency 
equations (which become power series in ^) and the 1 < p' < p corrections have been deter- 
mined. 

Eqs. (4.30/31) represent, of course, only a small subset of the full set of equations for the seven 
vertices (1.14). To establish the full set of self-consistency conditions by isolating the divergent parts 
of all DS interaction terms for these vertices represents, even at the lowest nontrivial (r = 1 and I = 1) 
level, a substantial research program, on which work is in progress. The example of (4.30/31) should 
however be sufficient to demonstrate that those divergent parts are considerably richer in structure 
than the perturbative ones, and that in this sense the evaluation with nonperturbatively modified 
diagram elements, rj^], squeezes more dynamical information from the DS dressing functionals. 



5 Special aspects of the quasiparticle subsequence 

Among the dynamical possibilities opened up by the extended iterative scheme, the "quasiparticle" 
subsequence - with r odd and only complex-conjugate singularities in the elementary two-point func- 
tions Dt and Sf - clearly is the most interesting one for QCD. Yet its complex singularities seem to 
raise special questions of interpretation and of possible violation of physical principles. While some 
of these were mentioned in a purely r = 1 context in the second of refs. [7], the crucial ones can be 
addressed adequately only with the perspective provided by the full sequential approximation. In the 
following we therefore discuss these conceptual questions one by one under appropriate catchwords. 

5.1 Short-lived elementary excitations 

The spacetime propagation characteristics implied by the two-point vertices (2.18) and (A. 25) of the 
quasiparticle subsequence are exhibited by continuing the corresponding propagators, D^' ' and Sp'°\ 
to Minkowskian {ku) space and Fourier-transforming to Minkowskian xm space (we recall that the 
order of these two steps is essential [30] if one is to arrive at the time-ordered Minkowskian two-point 
functions). We write results only for the scalar function Dt{x 2 m ), since the qualitative points will be 
the same for the two invariant functions of Sf(xm)- Working from (2.23), one has 
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where K\ is a modified Bessel function. By eq. (2.24), i.e. by virtue of asymptotic freedom, the 
small-distance behavior of this expression is that of a free propagator. At large timelike distances 
{x° M — > ±00 or x 2 M — > +00), assuming that s = refers to the pole pair closest to the origin of the k 2 
plane, the asymptotic behavior is 

Dp% 2 M ) -» Xr ' oA " 2 , e -WS£co6 L,o^+^o) , (5.2) 
with real parameters defined by 



n * (5-3) 

0r,2» 2 (crr,3) 4 = Xr,0 e = [ffr.O* 2 (^r,l ) 4 ■ 

In contrast to the "radiative" behavior of a real- mass propagator at timelike distances, this function 
decays exponentially, with an inverse lifetime, r _1 = 7 rj o, controlled by the imaginary parts of the 
leading /c 2 -plane pole pair. Note that the function decreases in both directions of the time axis. It 
describes an elementary excitation of the (transverse) gluon field, A%,(x)\Q >, that can exist only as 
part of a microscopically short-lived intermediate state. 

5.2 Fields without free-field limits 

The full Minkowski-space Green functions for N elementary QCD fields, denoted generically ipi with 
two-point functions Aj, of which at least one is a transverse-gluon field with propagator (5.1), have 
structure 

G N (x! . ..x N ) = < 0\T{(p!(xi) . . . A T (x k ) . . . (Pn(xn)} |0 > 

= D T (x k - xi)G N - 2 (. ■ •) + D T {x k - xi)D T {x m - x n )+ 

+ (permutations)] Gat_4(. . .) + (other disconnected terms) 

(5.4) 

+ / d A yi d 4 y N A 1 (x 1 -yi)... Dr{xk - Vk) ■ ■ ■ & N (x N - y N ) 
xf N (y 1 ...y k ...y N ), 

where Tjy is the connected and amputated function. For x° — > ±00 with the other x« fixed, the 
disconnected-gluon terms, if any, vanish exponentially by (5.2). The behavior of terms with only 
connected gluons, as exemplified by the fully connected last term of (5.4), is also dominated by 
relation (5.2) for the Dr(x k — y k ) factor, provided the y k integration is sufficiently convergent. The 
last point involves a subtlety: as discussed in sect. 3 and the appendix, the nonperturbative vertices 
from which Tjy is built have one real-axis pole (at k 2 M = -u^A 2 ) in the squared momentum of each 
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external gluon, and such poles would seem to cause insufficiently convergent, "radiative" behavior of 
the Tn in (5.4) with respect to (yk)M- However, as emphasized in connection with (4.8), the DS self- 
consistency conditions see to it that propagators acquire zeroes at the positions of vertex poles, so that 
the unamputated function, with propagators on all legs, retains only the complex momentum-space 
propagator poles, and therefore overall exponential-decay factors as in (5.2) for all transverse-gluon 
legs. The argument therefore goes through - the connected term, too, vanishes exponentially as 
x\ — > ±00. 

Now if all Green functions involving an At{x) field vanish in this way, then all matrix elements 
of At(x) between normalizable states decay exponentially as x° — > ±oo: the field has weak limits of 
zero at timelike infinity, in contrast to fields with stable-particle, real-axis propagator poles that are 
known to have nonvanishing weak limits representing free fields [31]. 

A different perspective on the same property results from looking at the quantity which "normally" 
- i.e. for fields with nonvanishing weak free-field limits - defines S-matrix elements: the r.h.s. of the 
reduction formula. This is obtained from the connected piece of (5.4) by going back to Icm space and 
conducting a residue search for all squared momenta approaching real mass-shell values: 

(S - H)(ki, . . . , k k , ■ ■ ■ , k N ) oc 

J™ 2 • • • J™ 2 • • • ,2 lim 2 { ~ ■ ■ ■ ( k l ~ t4) ■ ■ ■ fa* - tft) X (5-5) 

k x — k k -^n k k N ^n N 

x Ai(fc?) . . . D T (k 2 k ) . . . A N (k 2 N )T N (h, . . . k k , . . . k N )] } . 

In the "normal" case, the full propagators Aj have particle poles at mass shells kf = and the residue 
search locks in at those points to give a nonzero S-matrix element, proportional to the all-on-shell T. 
For a gluon with propagator (5.1), Dt(^) has only complex pole pairs staying at least a distance 
(Ima rj i)A 2 away from the real k\ axis, and by its numerator zero at k\ = i^A 2 cancels the real-axis 
pole of T. The residue search at real k%, no matter what the value of the real therefore gives zero. 
(The longitudinal gluons do retain their real-axis poles at k 2 = 0, as do ghosts in the present scheme, 
but these unphysical modes are of course excluded from the S matrix by appropriate projectors not 
explicitly indicated in (5.5).) One concludes that S-matrix elements involving at least one external 
gluon (or, by extension, quark) are, depending on semantic taste, either nonexistent (if one considers 
them as defined only for fields with nonzero free-field limits) or zero (if one applies definition (5.5) in 
an extended sense). The physical implication is the same - gluons and quarks are not asymptotically 
detectable reaction products. 

5.3 Violation of Causality? 

To avoid confusion here, we must recall that arguments deducing absence of complex singularities in 
amplitudes from microcausality refer to S-matrix elements (the propagators themselves, like all off- 
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shell Green functions, are always "acausal" by construction, i.e. nonvanishing at spacelike separations) 
and to those of the elementary fields (composite hadronic fields, being nonlocal, have no strictly 
microcausal commutators). The question then arises whether one-particle reducible contributions to 
S-matrix elements, in which one of the external sums-of- momenta flows through a single gluon or 
quark propagator, violate causality through their complex poles. (The simplest amplitudes of this 
type are the tree graphs of fig. 3). The answer is no, since these amplitudes also have external legs 
corresponding to the elementary (gluon and quark) fields, and therefore are zero (or must be declared 
nonexistent) as discussed in 5.2 above. In other words, the same set of propagators that seem to raise 
causality problems when occurring on internal lines, also ensure the vanishing of those amplitudes 
taken as S-matrix elements by being present on the external lines. 

5.4 Reflection Positivity? 

In Euclidean-coordinate (xe) space, the propagator Dp°\x E ) has an asymptotic form similar to (5.2) 
but with the roles of u> r fi and 7 r> o interchanged. At r = 1 and generally for low r in the quasiparticle 
sequence, such a propagator is not a purely positive function of x E , as required by reflection positivity. 
Its graph, shown schematically in fig. 4, has negative "overshoots" caused by the oscillations of K\ 
Bessel functions at complex arguments. This problem is potentially serious, but not inevitably so. 
It is quite possible, and indeed likely (given the fact that the DS equations are compatible with 
physical requirements like positivity) that it gets cured gradually, by the coherent admixture of faster 
oscillations from more distant complex-pole pairs, as the order r is increased. (This in fact is the 
normal situation in any discrete Fourier approximation of a positive function.) If this occurs, then 
the lack of reflection positivity at low r should simply be viewed as one of the natural errors expected 
in the low orders of any approximation scheme. It will basically be no more serious than e.g. the 
well-known violation of unitarity by the Born approximation, i.e. by low-order perturbation theory. 
However, it is true that such a conjecture can be checked only by an actual study of the r > 3 levels. 

If present, restoration-by-interference of reflection positivity would highlight a deeper point about 
the short-lived elementary excitations: the single complex-pole pair of the r = 1 propagator makes 
sense, strictly speaking, only as the lowest level of an improving sequence of approximations. While 
a stable-particle pole makes physical sense in isolation (describing then a free field), a quasiparticle 
pole pair can support itself only through the presence of interactions, which manifest themselves as 
dressing cuts. It is tempting, though at present speculative, to view this in parallel with the fact 
that a nonabelian gauge field cannot be defined without self-interaction. The formation of short-lived 
elementary excitations may be a dynamical specialty of nonabelian gauge theories. 
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5.5 Asymptotic Incompleteness 

The short-lived gluon or quark excitations as described by a two-point function in the " quasiqarticle" 
subsequence are distinct from ordinary decaying particles or resonances. The difference can be made 
precise by looking at the singularity structure of the propagator: for a resonance, this is character- 
ized by a physical, gauge-fixing independent cut on the real axis, whose discontinuity shows some 
local enhancement around the resonance energy. This indicates that the resonance will decay into 
stable, asymptotically detectable fragments, which together carry the open quantum numbers of that 
resonance. Therefore, it can in principle be recreated by macroscopically preparing beams of those 
fragments and colliding them. The propagator of the short-lived QCD excitation has no such decay 
cut. (Again, there is a subtlety - these propagators do have some real-axis cuts, arising from vir- 
tual conversion into longitudinal gluons and ghosts, the unphysical degrees of freedom, but these are 
gauge-fixing artefacts recognizable by their £ dependence.) The excitation therefore cannot decay into 
asymptotically detectable configurations carrying its open quantum numbers, and conversely cannot 
be prepared in isolation by colliding macroscopically preparable fragments. 

What can be produced from asymptotically controlled initial states are two-gluon and multi-gluon 
(or multi-quark) configurations compatible with the quantum numbers of the production channel, 
whose propagation in lowest order is described by product propagators of type 

{D(k 1 )D(k 2 )...D(k n )} a , (5.6) 

the notation indicating a coupling to the quantum numbers a. These propagators see to it that the 
configuration propagates only for a limited time of order A -1 before the dynamics of the theory forces 
it to pass its conserved quantities, and share of probability flux, back to some asymptotically accessible 
exit channels. 

The nature of the state space of such a system can at present be inferred only by educated guess. 
(Although the state space is in principle determined by the full set of Green functions, the prescription 
to construct the superficially convergent functions by skeleton expansions is far too inexplicit to 
actually provide a handle for an application of the reconstruction theorems.) Fig. 5 is a schematic 
representation of such a guess. It depicts the passing of unit total probability through the continuum 
space Ha of an ordinary scattering system, where it suffers some delay called a resonance, and through 
a subspace of exactly conserved quantum numbers, having an asymptotically accessible portion Ha, 
in QCD. In the latter case, probability flux can partially and temporarily - for times of order A -1 - 
leak out of Ha and populate a "closed", not directly accessible space, He- The presence of the latter 
implies that Ha does not, together with the bound-state space Hb, exhaust the full state space, i.e. 
that such a system lacks asymptotic completeness, a feature unfamiliar but by no means unacceptable 
physically. (Since temporal evolution in such a system maps a smaller onto a larger space and vice 
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versa, it may have to be described by operators that are isometric but not unitary, a feature not 
excluded by the requirement of probability conservation.) 

It is natural to inquire about the evolution of a (purely hypothetical, because of the impossibility 
of macroscopic preparation) single gluonic excitation, which would live in a subspace of conserved 
quantum numbers of the pure He type, without any Ha and Hb portions. This evolution would 
exhibit some analogy with that of a free particle in quantum mechanics, whose wave function, as long 
as no position measurement is performed, spreads indefinitely over all space. In a similar manner, the 
single short-lived excitation would convert into more and more complex multiexcitations, all short- 
lived, and thereby spread indefinitely over the whole He space. 

5.6 S-Matrix Unitarity 

While evolution projected onto the asymptotically accessible subspace Ha is not unitary at small times, 
due to the partial leakout of probability into the He space, the fact that the leaking is dynamically 
limited to times of order A -1 guarantees that between asymptotic times, t = -co and t' = +00, 
probability conservation does hold on that subspace alone, since all probability flux fed into the 
system through scattering eventually gets pushed back into asymptotically accessible exit channels. 
In other words, one expects the S matrix to remain strictly unitary. 

It is instructive to check this expectation out in a more formal way by asking whether the relation 



is compatible with the presence of the "closed" subspace He- If at least one of the two states \a >, \ b > 
is in that subspace, then by the discussion around eq. (5.5), the l.h.s. and at least one of the factors 
on the r.h.s. must be set equal to zero, and the relation is trivially fulfilled. If both \a > and \b > 
are in the asymptotically accessible space Ha, then states \n > in He, again by the arguments of 
sect. 5.2, make no contribution to the sum on the r.h.s. But then the imaginary part on the l.h.s., 
which is equal to the discontinuity across the timelike real axis of the relevant energy-squared variable 
for the a — > b process, does not receive contributions from such states either: intermediate states in 
He are described by propagators of type (5.6), and integration over them produces pairs of cuts at 
complex-conjugate locations, which together always leave the amplitude continuous and real along the 
real s axis (even though some of them do start at points on that axis), i. e. which do not produce 
physical absorptive parts. ( This can be checked out in detail for e. g. the one-loop terms of fig. 2 
with the r = 1 propagators ). One concludes that it is indeed consistent for eq. (5.7) to hold on the 
asymptotically accessible space alone. 

The discussion of this section may be summarized by stating that, perhaps contrary to first ap- 
pearances, a solution {r[ r ' p l|r odd} of the quasiparticle type with complex-conjugate propagator sin- 




(5.7) 



n 
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gularities does not necessarily imply violation of physical principles beyond the natural errors expected 
in the lower orders of any systematic approximation scheme. Such a solution is in fact well suited for 
a description of the elementary QCD excitations where they come closest to being "seen", namely at 
the origins of jet events. 
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Appendix 



A Other Superficially Divergent Vertices 
A.l Four- Vector Vertex 

This is the highest of the superficially divergent vertices (1.14). Its complexity with respect to color 
and Lorentz structure exceeds by an order of magnitude that of all the other T sa n v combined, and is 
the main source of technical difficulty in the p = self-consistency problem. The number of color- 
basis tensors - each a fourth-rank object over the adjoint representation of SU(Nc) - in the 
decomposition 

L 4V,abcd — ^abcd 1 4V(i) ( A - L ) 
i 

is eight for Nc = 3 and nine for Nc > 4. For brevity, we refer the reader to ref. [11] for the general 
case and write an example of a suitable color basis only for the case of immediate interest in QCD, 
viz. Nc = 3: 

Cibcd = dab5 c d , Cabcd = &achd , C^L = &ad&bc , 

C 'abcd = ■) ? abn j "cdn i ^ a bcd = f acn fd^ n ~ f adn f ben i (A- 2) 

C a bcd = dabnfcdn , C^feL = dacnfdbn , CabL = ^adnfbcn ■ 

The combination + is not listed, as it equals —C^ by virtue of the Jacobi identity. In the 
many situations where minimality of the basis is not crucial but manifest Bose and crossing symmetry 
are, we shall nevertheless use the nine-dimensional set including and separately. 

For Lorentz structure, we again concentrate on the totally transverse pieces defined in analogy 
with (3.2), which have tensor decompositions 

rgg v (p! . . . P4 ) = t«'% Pl )t x '\p 2 )t»'^p 3 )t^(p,) 

(A.3) 

x E E L liT)(Pi ■■■Pa) G {i)kj {p\...p\\ s, t, u) , 

k=o j K ,J ' 

with J2Pi = again understood. Linearly independent Lorentz tensors constructed from the 

Euclidean metric 5 alS and three independent combinations (denoted generically ^1,2,3) of the pi, come 
with three mass dimensions, 2k = 0, 2, and 4. They include three dimensionless (k = 0) tensors, 

LfSg = 5 kX 6^ , LjJ^ = 8 KIX 5 Xu , Lffifi = 5 KI/ 5^ , (A.4) 

plus 54 independent 2k = 2 tensors of type S^k^k^ (to, n € {1, 2, 3}), of which 24 contribute to 
the totally transverse vertex, plus 81 independent 2k = 4 tensors L( 2 j) of type k^kfk^k^ (k, l,m,n € 
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{1, 2, 3}), of which only 16 contribute to the 4T vertex. The perturbative zeroth-order vertex contains 
only the k = Lorentz tensors (A. 4) and has the expressions 



r (0)pert =c(4) 



L 



(0,2) 



L 



(0,3) 



+ C( 5 ) 



L 



(0,3) 



L 



(0,1) 



+ C( 6 ) 



L 



(0,1) 



(0,2) 



=cw [1 ( 



(0,2) 



(0,3) 



+ 



i (c( 5 ) - c( 6 )) 



(A.5) 



-2L 



(0,1) 



+ ^(0,2) + ^(0,3) 



in the linearly dependent but crossing symmetric and in the minimal color basis (A. 2), respectively. We 
do not write full enumerations of the k > objects here, which will be of interest only in connection 
with detailed loop computations. For a first exploration of the solution described here, it will be 
reasonable in any case to start by looking for approximate self-consistency with a strongly restricted 
tensor structure. 

One would expect to treat the invariant functions G^^j associated with these tensors as functions 
of six independent Lorentz-invariant variables, an example being the six scalar products formed from 
the three conserved total four- momenta in the three crossed two-gluon channels, 



P 
R 

Q 



Pi +P2 
Pi +P3 
Pi +PA 



~(P3 + Pi) , 
~(P2 +Pi) , 
~{P2 +P3) • 



(A.6) 



Again, however, it turns out that one has no freedom here: dynamical consistency in the three- 
gluon DS equation connecting T^y to T^y, together with crossing symmetry, dictates that rational 
approximation of G's be performed with respect to the seven variables indicated in (A. 3), where 



s = P 2 



u = R 



t = Q 2 , 



and with the relation expressing their linear dependence, 

4 

s + u + 1 = p. 



2 

i i 



(A.7) 



(A.8) 



i=i 



being carried as a subsidiary condition. This would seem to considerably complicate the writing and 
proper restriction of FDRA's were it not for two simplifying features which we anticipate: the various 
additive terms of (A.1/A.3) can be grouped as either crossing triplets or crossing singlets, and if one 

\r 01 

restricts attention to T^ , the nonperturbatively modified vertex of zeroth perturbative order, then 
self-consistency is possible for a simplified structure, where in each member of a crossing triplet the 
invariant-function denominator contains only one of the Mandelstam variables s, t, u at a time, while 
for a crossing singlet it does not depend on them at all. Moreover, pole positions in FDRA's with 
respect to the four pf and with respect to the variables (A.7), which a priori could be chosen differently, 
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are forced by DS consistency to be in fact the same. Thus the invariant function for the s-channel 
member of a crossing triplet would possess r-th degree FDRA's of the form 



r [r,0] 



^4v(i)fcj(Pi---^; s ;5(«-*)) 



jft + < 2 A 2 )] { n [ft {pf + < 2 .a 2 )] } 

and the other two members of the triplet would follow from this by the two crossing operations 



(A.9) 



(b,X,p 2 )< — ►(c,/i,p 3 ); P^R,s^u; 
(b,X,p 2 )< — > (d, f,p4); P^Q,s^t. 



(A.10) 



The mass dimension —2k, as well as the various restrictions and Bose-symmetry properties, are built 
into the numerator polynomial, 



N iv(i)k,j (pi ■■■Pl,s;±{u- tfj = (S u + S i5 + S i6 ) 5 k0 (spIpIpIpIY 



^ n (i)k,j(r) 



jQ (p2) m ! s n i (Hzz*)" -2 ^2^5r-fc-(mi+...+ri2) _ 

:=i 



(A.ll) 



mi...«2>0 

Condition (2.3) allows nonzero coefficients C mi ... n2 only for 

mi + 7TT-2 + 777-3 + 7774 + 7li + 772 < 5r — — 1 



(A.12) 



The additional restrictions arising from condition (2.5) are now more involved: when T^y appears 
in a 1PI diagram, up to two of its four legs may be external. To avoid the occurrence of ultraviolet 
divergences stronger than the perturbative ones, no term of the vertex should then behave worse than 
a constant at large loop momenta when (i) any three of its momenta pi are running in loops and 
(ii) any two of its momenta are running in a loop. Considering only the case where the 2k momenta 
supplied (for k > 0) by the tensor are all running along in those loops, we find that (i) imposes 

the four restrictions 



^(1 - 6 H )mi 



+ ("I + n 2 ) < 4r - k (i = l,2,3,4), 



(A.13) 



whereas (ii) implies, for the function (A.ll), the six restrictions 

(777j + rrij) + (n\ + 712) < 3r — k 

= (1,3), (1,4), (2, 3), (2, 4); 

2(tt7i + 777 2 ) + 77 2 < 4r - 2k ; 
2(7773 + "74) + "2 < 4r - 2k . 



(A.14) 



(A.15) 
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In writing (A. 15) we have used the fact that 



u - t = (p 2 -pi) • (p 4 -p 3 ) . (A.16) 

In other cases where some of the "tensorial" momenta stay out of the loops, these restrictions may be 
relaxed slightly, but the additional freedom gained in this way is minor. 

For the self-consistency problem, the relevant structural property emerging from the above is that 
with respect to the variables (A. 7), the p = vertex may be assumed to have a simplified partial- 
fraction decomposition which is the analog of (3.10), 



71=1 



p(r) ( A 2 , 
E n(s) [7T^ ' 



(A.17) 



'"(«) I u+u" A 2 / ^ "(*) I A 2 



Here E n r u \, E n u\ are the crossing partners of E n t s \ in the sense of (A. 10), and each of the E n tensors 
in turn has a partial-fraction decomposition with respect to the pf variables whose form is exemplified 
by 



l ar fir / ^ \ a / ^ y 

E n(s) = E EE E n(s)ar o „ A2 2 „ A2 • (A. 18) 



a,/3=0 o-=a T= (3 

Here the E n ^ UT functions still have dependence on p%p\. 
A. 2 Ghost vertices 

The well-known fact that diagrams, whether perturbative or dressed, with external G and G lines 
factor out one power of an external momentum per GG pair, causes the Dyson-Schwinger equation 
(1.13) for the inverse ghost propagator to assume the special form 

5 ab r G e(p 2 ) = -<W 2 jl+ l[D,D,T Gvls ] (p 2 )| , (A.19) 

where the dimensionless, logarithmically divergent loop integral I(p 2 ) involves the Tgyg three-point 
vertex. Rational approximants of the odd-r sequence for the ghost propagator D are therefore of a 
type we omitted as exotic when discussing Dt — they generically have two different real-axis poles, 
one massless and one massive, which for an unphysical excitation is not impossible. For example, the 
counterpart of (2.15) reads, 

D^\ P \^)= 2 P ^y A ly (A.20) 

This includes the special cases = ui,i, where the r = 1 function reduces to an r = function with 
one massless real-axis pole, v\ : 2 = 0, where the same happens but with a mass (generally gauge-fixing 
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dependent) for the FP ghost, and v\ : \ = 0, where the ghost develops a second-order pole at p 2 = 0, as 
it does under the Gribov-Zwanziger mechanism [14,12]. For a generic gauge fixing, however, (A. 20) 
and its higher-r extensions are likely to remain the relevant forms. 

In the ^q V q three-point vertex, FDRA's for the invariant functions F^(c = f or d) defined by 
the tensor decomposition 



r a£(Pi> fc > -Pi)gvg = fo.bc PiF(f)o + k"F {f)1 



+ d, 



abc 



vtF [d)0 + k?F m 



(A.21) 



(where now p\ + k — p 2 = 0), are somewhat more complicated than (3.6) because of the absence of 
symmetry properties. In particular, denominator zeroes in the ghost-line variables p\ and p 2 , which 
we denote as - f^ 2s A 2 , may differ from those in the gluon variable k 2 even at r = 1. The general 
structure is now 



F, 



y,o] 



(c)i 



n (k*+K,2s^) 

8=1 



n (pi+K^) 

3 = 1 



s=l 

(c = f ot d, i = or 1, r = 1, 3, 5 . . .) , 
where numerator polynomials must respect the perturbative limit 

^(0)pert 



(A.22) 



(c)i 



as well as restrictions analogous to the k = version of (3.8/3.9), but no symmetry restrictions. The 
explicit form of the i ? (/)o function for r = 1 will again be given for illustration: 

F lf$ (Pi fc2 '#D = 1 + Vi%i + V3% 2 + 

+ (w + y7%i) n Pl n P2 + n fc [ V2 + (y 4 + ys $) n P2 (A.23) 
+ (2/8 + V9j?) n P i + 2/ion pl n P2 

Here II P = A 2 / (p 2 + f 2 A 2 ) denotes a rational building block analogous to the II of eq. (3.11). The 
important point for the DS self-consistency problem is again the existence of a "regular-plus-singular" 
decomposition, analogous to (3.10), in any one variable; e.g., 

[r,0] 



T GVG(c) iPl ' k >~ P2) 



- B» c)0r (p 2 ,p 2 ;k 2 ) 
+ E Pi) 



s=l 



A 2 



(A.24) 



In refs [7], using heuristically what we would now call the r = 1 degree of rational approxima- 
tion, treatment of the two ghost vertices was subject to the prejudice, taken over from older work 
of Eichten and Feinberg [9], that vertices do not develop nonperturbative corrections with respect to 
the momenta of unphysical external lines. This, in particular, kept the ghost propagator iM 1 ' ! in its 
perturbative form, 1/p 2 . From the more systematic viewpoint advanced here, this prejudice is seen 
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to be unjustified. Even if DS dynamics wouid make the case -17,2 = iV,i prevail for all r, the function 
iM 7 *' !, which in contrast to the longitudinal gluon propagator Dl is not protected by an ST identity, 
would pick up pole-zero pairs at r > 3 to represent its dressing cuts, and these in turn would demand 
p\ and p\ denominator factors in (A. 22). This is in accord with the fact, emphasized by Lavelle and 
Schaden [24] in an OPE context, that nonperturbative ghost-antighost vacuum condensates of zeroth 
perturbative order do exist. 



A. 3 Fermion propagators. 



We use Euclidean 7 matrices obeying {7^,7^} = — 2<5 /w/ , so that jyp = —p 2 . The flavor-F, inverse 
propagator of zeroth perturbative order, analogous to (2.18), has odd-r approximants 



,(F)\ 



(k {1 I 



(r+l)/2 

+ £ 

r,2 s=2 



(k (F) V (k (F) V 

\ K r,2s+l J \ K r,2s+l J 

7742 + 774? 



(A.25) 



Here we encounter a complicating feature: the nonperturbative mass scales n r ,i cannot, except in 
special cases, be written simply as numerical multiples of A, since there now exist additional RG- 
invariant mass scales. These are the extraneous, or Lagrangian, mass scales rhp, one for each flavor 
F, that are connected to renormalized quark masses mjr(y) through 



(m F ) R = [m F (u)] R ex.p < 



dg' 



7m ( 



/) 



R 



(A.26) 



in the scheme indicated by R. The rational character of approximants with respect to A then allows 
fermionic mass scales, such as the k's and of (A.25), to be polynomials in A and rhf. The only 
restrictions come from condition (2.3), which demands that the bare vertex, 



be approached in the "perturbative" limit, A — > 0. In general, we then have 

(F) . 

rrif + A , 



JF) 
K r,l 



(F) '(F) * (F) . 

<,2s = w r,2s m f + w r,2s A . 



(k (F) V 

\ K r,2s+l J 



A 



f (F) ( 



(A.27) 

(A.28) 

(A.29) 
(A.30) 



The primed coefficients, w', are absent in the strict chiral limit, rhp = 0. In the opposite limit, 
t%f » A, the pole terms of (A. 25) are of order A/rhf relative to the leading rhf term of (A. 27). To 



41 



illustrate (A. 25), we write the analogs of (2.13) and (2.15) for r = 1 (omitting now F on mass scales): 



T 



[1.0] 
FF 



2 

1,3 



t 2 



1.2 



+ n 



*l,2*i, 3 
K l,l + „2_i_^.2 



(A.31) 



-,[1,0] _ i±JS£ 

'F - 

-^[p 2 + (ki,2(ki, + +ki,_)-|ki, ± | 2 )]+H[(ki, + +ki,_-ki,2)p 2 +ki,2|ki,±| 2 ] 

_ (P 2 +K 2 + )( P 2 +K 2 _) 



(A.32) 



where II denotes a unit Dirac matrix, and where 



«i,± 



5 K 



1,2 



K 



1,3 



"1,3 



1 , 

2Kl 



K l,2, 



(A.33) 



A complex-conjugate pair of poles, and therefore a short-lived quark-like excitation, is present if 



«i,i + K i,2 7^ and /cf 3 > 



2 («i,i-«i,2) 



(A.34) 



Again, for the interpretation to hold at r > 3 requires a definite pattern of one nearest-to-origin pole 
pair, plus well separated zero-pole strings identifiable as complex cuts, to emerge as r is increased. 

It is interesting to note the tight correlation present in (A. 31) between spontaneous chiral-symmetry 
breaking and the emergence of short-lived quarklike excitations. In the chiral limit rhf = 0, where the 
Ki : i scales are purely spontaneous (multiples of A), chiral-symmetry breaking can be avoided only by 
having the invariant function associated with the Dirac 11 in (A. 31) vanish identically, which implies 
either = Kip = 0, or ki ; i = ftf 3 = 0. In both cases, the pole positions (A.33) turn purely real. In 
other words, whenever short-lived quarks are formed dynamically in the rhf = case, chiral symmetry 
must be broken. 

The second line of (A. 31) moreover illustrates a technical but important point: residues of the two 
invariant functions of Sp' ^ 1 at the commmon p 2 = —k\ 2 P°l e cannot be chosen independently but 
must be related in such a way as to combine into a single pole in the matrix-valued variable fi, 

-$+Ki2 1 



p 2 + k\ . 



■fi + Kl, 



(A.35) 



The reason is that by using unrelated residues one obtains, in general, a propagator Sf with three 
poles, and therefore ends up in the subsequence describing stable quark particles that we do not wish 
to study here. 
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A. 4 Fermion-vector three-point vertices. 

Color structure is strictly identical to that of the bare vertices: 

r K«'^-(5 A *L^>'*'- w) (A - 36) 

(c, c', triplet-color indices; a, (3, Dirac-Spinor indices, k, gluon four-momentum, and p\ + k — p 2 = 0). 
Lorentz (matrix- vector) structure of the amplitudes is then the same as for the QED vertex [32]. 
We therefore mention only that a tensor decomposition suitable for our purpose, 

12 

Kp(pi, k, -P2) = £Wi (pi k 2 ,p 2 ) (V?) ap , (A.37) 
1=1 

is possible in terms of the twelve matrix-valued vectors, 

V 1 =-f, V 2 =£ i r, V 3 = -ffi 2 , V A =hrf> 2 , (A.38) 

V 5 = r^, V 6 =^, V 7 = r^ 2 , V a =fi 1 r»j 2 , (A.39) 

V9 = fc"n, F 10 =A^, F n = ^ 2 , V 12 =^k^ 2 , (A.40) 

where 

r" = ^(p2+Pir, ^ = (P2-Pi)". (A.41) 

(The vertices for transverse gluons, T FT p, clearly contain only the first eight, eqs. (A. 38/39), of these 
matrices). In this basis, the restrictions imposed by C invariance, which leave nine of the twelve scalar 
functions W% independent, take rather simple forms. In particular, 

Wi (p\, k 2 ,p 2 ) = W i+1 (pi k 2 ,p 2 ) (i = 2 or 6) , (A.42) 

W i {p 2 ,k 2 ,p 2 )=W i {p 2 ,k 2 ,p 2 ) (i = l,4,5,8), (A.43) 

so that the transverse part has six independent invariant functions, with four of these being symmetric 
under interchange of the fermion and anti-fermion variables. 

Construction of the FDRA sequence is greatly simplified by anticipating that DS self-consistency 
will dictate, for the fermionic variables, the correlation between pole and matrix structures indicated in 
(A. 35), leaving only " matrix- valued poles", with the ordering of matrices adopted in (A.38 - 40) always 
observed. The resulting transverse, nonperturbative vertex of zeroth perturbative order, obeying 
conditions (2.3/2.4), will be written in the form of a partial-fraction decomposition: 

r^pCpx, k, - P2 )} M1 = jc^fo, -p 2 ) + £ C&.CP!, - P2 ) ( fc2+ ^ A2 ) } > (A-44) 
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where 



Co!r(Pl'-P2) = 



s=l 



r,2s 



-i- v a r (o ^ , JO 



tiaras' 



(A.45) 



+ 



terms with (k 



2\n 



1< n < 



r-l 



^oV + ^ ( o r) (5f)+<i r) (x7^ + 7^) + 

+ £ ^ [(41 + *2&£) 7" + 41^ + 4t£ 

'(r) r^ol , ^ \( (r) , '(r) p? \ „ . (r) j/i u 
s=l 

A t h J2 A - 



(A.46) 



=1 j/i+*r, 



,7 M + 2 ; 



(r) r M 



A 



i/2 r 2s' 



Again, preservation of perturbative renormalizability turns out to be a requirement slightly stronger 
than (2.3/2.4), imposing extra restrictions in the case where p\,P2 run to infinity in a loop while k 
stays outside the loop and fixed. For the vertex to behave no worse than a constant at large loop 
momenta in this situation requires the primed coefficients in (A.46), or sums of these, to vanish: 



Jir) _ »(r) _ 



%0 



'(r) _ 



Hfl 



E 



'i,3,s — 



(1 < t < r, all r) ; 

(1 < t < r, all r) . 



(A.47) 
(A.48) 



s=l s=l 

In the remaining terms, mass scales have for simplicity been written in a form appropriate for the 
chiral limit, rhf = 0. At rhf ^ 0, those mass scales in (A. 45/46) that are not forced by condition (2.3) 
to be pure multiples of A should be generalized in an obvious way, following the pattern of (A. 28/29). 
Formally, this can be accounted for by allowing for a dependence on the dimensionless ratio rhf /A in 
the relevant dimensionless z coefficients. In fact, in sufficiently high loop orders, all nonperturbative 
coefficients, including those of the gluon-ghost sector, will eventually depend on all rhp/A ratios. 
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Figure Captions 



Fig. 1 Distribution of poles (crosses) and zeroes (circles) for nonperturbatively modified propagation 
functions Dp°\k 2 ) in (A) "particle" subsequence, (B) " quasiparticle" subsequence. 

Fig. 2 Diagrammatic form of Dyson- Schwinger equation for two-point vertex (negative-inverse prop- 
agator) of gluon field, showing one-DS-loop terms (A . . . D) and two-DS-loops terms (E&.F) of 
dressing functional. 

Fig. 3 One-particle-reducible contributions to four-point amplitudes of elementary QCD fields that 
seem to produce causality violation in S-matrix elements for the "quasiparticle" subsequence. 

Fig. 4 Qualitative behavior of propagation function D in Euclidean coordinate space for free massive 
propagator (narrow line) and a "quasiparticle" propagator of type iM 1 ' ! (broad line) 

Fig. 5 Schematic view of probability flow through state space of (I) an ordinary scattering system 
and (II) an asymptotically incomplete system with a closed subspace He- 
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(WWW) = - ( wwvwwv ) + 
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